Fungal Genetics and Biology

Fungal diseases cause devastating crop infestations and millions of deaths annually. To improve disease treatment development, we must understand the biological processes these treatments target and the fungal cellular responses they elicit. Current methods to identify genetic responses of fungi to antifungal drugs and fungicides typically involve laborious molecular genetics experiments or transcriptomics. We utilized high-throughput screens in the basidiomycete yeast Rhodotorula (Rhodosporidium) toruloides to rapidly identify genes associated with the fungal response to a copper-zinc fungicide developed by VM Agritech and compared identified genes to those transcriptionally activated in response to the fungicide in R. toruloides and the filamentous ascomycete Neurospora crassa. These 'omics technologies combined with molecular genetics and biochemical analyses showed copper and zinc synergize to inhibit fungal growth due, at least in part, to an increase in cellular zinc concentrations when copper is present. Intriguingly, although we identified genes involved in the copper-zinc fungicide response using both high-throughput screens and transcriptomics, there was little overlap in genes identified through the two methods. Our findings demonstrate the power of combining functional genomics and transcriptomics approaches in identifying a more comprehensive map of the genetic mechanisms necessary for responding to fungicides.

Introduction. Histoplasma capsulatum is an important fungal pathogen of humans, especially in immunocompromised patients, and the World Health Organisation has listed Histoplasma as a High Priority fungal pathogen. The taxonomy of the genus Histoplasma is still under development, but its only species H. capsulatum has at least of seven lineages, and a range of sublineages. Gap Statement. Allelic multilocus sequence typing (MLST) has been very successfully applied for typing and epidemiology of pathogenic bacteria, but its adoption for fungal pathogens is more limited. Genome sequence-based MLST approaches promise much a higher resolution for comparative phylogenetics, but there are only a small number of genome assemblies available for H. capsulatum. Aim. To develop and validate genome-based MLST typing for H. capsulatum. Methodology. A total of 406 genome assemblies were generated from publicly available Illumina short-read sequencing samples. Selected genome assemblies were subsequently used to generate MLST schemes with chewBBACA, based on whole genome sequences (wgMLST), NCBI-obtained annotations (annMLST) and annotations obtained with Funannotate (fanMLST), and used for lineage identification, and comparative analyses of population structure. Results. Genome assemblies were successfully generated for 406 H. capsulatum samples, with BUSCO completeness scores of over 98%. A subset of 60 genome assemblies was used to generate and validate MLST schemes based on the 60 genome assemblies, 5 genome annotations obtained from Genbank, and 13 genome annotations generated using Funannotate, resulting in schemes with > 8,000 markers. Each of the MLST schemes allowed reproducible separation of all 7 lineages and 6 H. capsulatum Suramericanum sublineages, without the need for high-performance clusters or powerful computing equipment. Comparison of the different schemes using tanglegrams and scatter plots showed good reproducibility (r2 values of 0.8-1.0). Comparison of total sizes of genome assemblies showed significant lineage-dependent size differences, with the Suramericanum, LAmB, India, Capsulatum and Mississippiense lineages averaging a total genome size of 30-33 million bp, and the Africa and Ohiense lineages averaging at 39 million bp. Conclusion. Genome sequence- and especially genome annotation-based MLST analyses are a promising approach that can be adapted as a rapid, high-throughput method for genomic clustering and epidemiological tracing of pathogenic fungi.

Abstract Tricyclic coumarins (TCs) exhibit defensive properties against herbivorous insects and possess valuable pharmacological properties for humans. They are mainly distributed in the Apiaceae, Rutaceae, Moraceae, and Fabaceae families. The chemical structural diversity of TCs in these families primarily arises from differences in skeletal variety and attachment position of the third ring, as observed in unsaturated furanocoumarins (FCs) and saturated dihydrofuranocoumarins (DHFCs). In TC biosynthesis, the FC and DHFC pathways diverge through distinct oxidation patterns. The product specificity of the reported cytochrome P450 monooxygenase at the branching point suggests that these TC groups are produced in nearly identical ratios; however, some species of the genus Angelica (Apiaceae) accumulate DHFCs at higher levels than FCs. In this study, we investigated the enzymatic mechanism responsible for forming these two TC groups using A. keiskei, a Japanese native species used as a medicinal and edible herb. Organ-specific transcriptome analysis highlighted a P450 member CYP71AJ62 that exhibited a close correlation with TC accumulation across organs. Enzymatic characterization demonstrated that CYP71AJ62 predominantly forms a hydroxylated DHFC, vaginidiol, which may promote DHFC over FC skeletal formation in TC biosynthesis. Phylogenetic analysis suggested that CYP71AJ62 has been derived from a gene responsible for FC formation by functional shift without gene duplication. Subsequent three-dimensional structural modeling and site-directed mutagenesis revealed that two amino acid residues on substrate recognition sites cooperatively determine the product specificity of this P450 enzyme. This study provides new insights into the molecular evolutionary process underlying ring-structure diversification and, more broadly, the adaptation of plants to their environments.

Macaranga tanarius (Euphorbiaceae), a myrmecophilic plant grown in subtropical countries, was identified as a source of Okinawan propolis, a honeybee product, which is glandular trichomes developed on the surface of the fruits. This tissue contains a high amount of characteristic geranylated flavonoids known as nymphaeols (nymphaeol A - C and isonymphaeol B), which exhibit strong antioxidant activity. In this study, we characterized the structure of glandular trichomes of M. tanarius, showing that they are sac-like structures containing drupelet-like aggregates of secretary cells, and a subcuticular cavity was also observed. From a cDNA library, we identified two prenyltransferase (PT) genes involved in biosynthesis of nymphaeols, and a PT gene grouped in a primary metabolic PT. Our biochemical experiments demonstrated that one of the former secondary metabolic PTs (MtPT1) is a B-ring-specific geranyltransferase for eriodictyol to yield nymphaeols B and isonymphaeol B. MtPT1 was unusual for a plant PT, an enzyme that yielded multiple reaction products. The other PT (MtPT3) transferred a geranyl moiety to the A-ring of eriodictyol to form nymphaeol A. Although dimethylallyl diphosphate (DMAPP) was not recognized as its prenyl donor substrate of eriodictyol, giving a B-ring geranylated eriodyctiol as the prenyl acceptor substrate, this enzyme used DMAPP to attach the prenyl moiety to the A-ring and produced nymphaeol C. These data suggest that M. tanarius has evolved its PT function to produce a diverse range of prenylated flavonoids using a limited set of genes.

Abstract: Picrotoxane-type sesquiterpene alkaloids, including dendrobine, are characteristic constituents of Dendrobium nobile, but the enzymes that produce their proposed skeletal precursors remain in completely characterized. Here, we combined tissue metabolite profiling with genome-guided screening and biochemical characterization of terpene synthases (TPSs). GC-MS analysis revealed variation in sesquiterpenoid profiles and dendrobine signals among tissues and developmental stages. Manual curation of the chromosome-level genome yielded 16 full-length TPS-a candidates. Functional screening in an FPP-supplying Escherichia coli system showed that three closely related genes on chromosome 1, DnTPS1~DnTPS3, produced copaborneol as their dominant product. The principal product of DnTPS1 was purified and identified by NMR spectroscopy, and a purified-enzyme assay confirmed its direct conversion of FPP to copaborneol. Fluorescence imaging in rice protoplasts and Nicotiana benthamiana leaves indicated predominantly cytosolic localization of DnTPS1. Guided by sequence comparisons and structural modeling, we tested 12 DnTPS1 variants and found that substitutions at W377 and I295 altered the relative distribution of sesquiterpene products. Phylogenetic and regional gene-order analyses placed DnTPS1~DnTPS3 in the same clade and identified related TPS candidates within corresponding genomic regions of six Dendrobium species. These findings establish that three D. nobile TPSs can form copaborneol and identify DnTPS1 residues that influence product selectivity, providing a basis for testing the proposed connection to dendrobine biosynthesis.

Pectobacterium, a member of the soft rot Pectobacteriaceae (SRP), can cause highly destructive soft rot and blackleg diseases on many important crop and other host plants, particularly potato (Solanum tuberosum). Pectobacterium and its sister genus Dickeya have common virulence strategies, yet a highly diverse "open" pangenome containing many unique and hypothetical proteins. Here, we used randomly-barcoded transposon insertion-site sequencing (RB-TnSeq) on six Pectobacterium species to identify genes involved in bacterial survival and growth in potato tubers, as well as in minimal media conditions containing different single carbon sources. The strains evaluated varied broadly in genomic content as well as host plant and year of isolation. We focused on identifying conditionally essential genes across all strains tested, as well as genes that were important for growth in one or few of the tested strains. These results highlight the importance of many common genes involved in metabolism and maintaining cellular homeostasis, as well as shared and rare virulence traits including polygalacturonate catabolism, genes involved in polysaccharide biosynthesis, immunity proteins, and (for P. brasiliense) self-tolerance of the antibiotic carbapenem. This large dataset of conditionally essential genes, generated across Pectobacterium species during soft rot infection of the agriculturally important plant host potato, can serve as a rich dataset for further hypothesis generation and identification of important and/or novel in planta fitness factors.

The Fusarium oxysporum species complex (FOSC) comprises numerous plant pathogens infecting crops such as tomato, lettuce, melon, and banana. Despite its economic and ecological importance and decades of research, the species names linked to different phylogenetic lineages within this species complex remain largely controversial. Earlier taxonomic frameworks relied on a limited number of marker genes, variable taxon sampling and inconsistent species definitions, leading to a fragmented and unstable taxonomy. To address these challenges, we here analysed nearly 2,500 conserved marker genes from 660 genomes, including genome assemblies from 22 ex-type strains that were sequenced using long-read sequencing technology, representing nearly half of the 51 currently described species within FOSC. Maximum likelihood- and coalescence-based analyses of the full phylogenomic data matrix as well as subsets yielded a robust phylogeny of FOSC, which was furthermore corroborated by whole genome comparisons. Many previously described species are redundant, while two novel lineages remain undescribed. Direct comparisons to traditional marker gene-based approaches demonstrate their inadequacy for species delimitation in FOSC. We show that our phylogenomic framework can guide the detection of novel marker genes able to distinguish the major lineages within FOSC, which will be essential for taxon delineation and classification in future studies. Thus, our phylogenomic analyses provide a stable cornerstone for revising FOSC taxonomy and offers a critical resource for the broader research community.

Background: Plant defensins are small cysteine-rich antimicrobial peptides. In Arabidopsis thaliana, the inducible Class I gene AtPDF1.2a is widely used as a transcriptional marker of jasmonate/ethylene immunity, whereas the constitutively expressed Class II members have received less functional attention. AtPDF2.2 was compared with AtPDF1.2a during infection of Arabidopsis by the vascular fungal pathogen Verticillium longisporum, the necrotrophic fungus Sclerotinia sclerotiorum and the bacterium Pseudomonas syringae pv. tomato DC3000. Results: Profiling of the AtPDF family showed that AtPDF2.2, together with the Class II genes AtPDF2.3 and AtPDF2.5, was significantly reduced during V. longisporum infection, whereas AtPDF1.2a was induced. A promoter-reporter line showed a corresponding decrease in AtPDF2.2 promoter activity during V. longisporum infection, strongest in infected roots. KD-pdf2.2 plants contained more V. longisporum DNA, whereas OE-PDF2.2 plants contained less. The same reciprocal disease pattern was observed during S. sclerotiorum and Pst DC3000 infection. Histochemical Prom::AtPDF2.2::GUS reporter signal decreased during S. sclerotiorum infection but increased following Pst DC3000 infection. Basal AtPDF2.2 expression was lower in JA and ET-signaling mutants and higher in the SA-deficient sid2 background, while the analyzed promoter lacked the GCC-box that mediates inducible AtPDF1.2a regulation. AtPDF1.2a overexpression reduced symptoms across the three pathosystems, but loss of AtPDF1.2a increased disease only during V. longisporum infection under the conditions tested and did not alter Pst DC3000 DNA abundance. AtPDF2.2 overexpression was not associated with reduced vegetative growth under the conditions tested, whereas AtPDF1.2a overexpression reduced rosette size and leaf number. Conclusions: AtPDF2.2 is a constitutively expressed Class II defensin whose expression level is associated with disease outcome across two fungal pathogens and a bacterium. Its promoter shows tissue-dependent regulation during V. longisporum infection, while histochemical reporter signal decreases during S. sclerotiorum infection and increases during Pst DC3000 infection. This regulatory behavior differs from that of the inducible reference gene AtPDF1.2a. The reciprocal disease phenotypes of the AtPDF2.2 transgenic backgrounds support a protective association that may involve direct antimicrobial activity, modulation of host responses or both. The absence of reduced vegetative growth in OE-PDF2.2 under the conditions tested supports further evaluation of PDF2-type genes in Brassica.

Extremophilic fungi are a relatively unexplored source of stress-resilient microbes for sustainable agriculture. Here, we describe Talaromyces pseudofuniculosus FAR, isolated from the highly acidic site of Piccolo Inferno in southern Italy and identified by ITS and {beta}-tubulin (BenA) sequencing. Despite its origin from an extreme environment, T. pseudofuniculosus FRA remained fully functional under non-stress conditions and displayed tolerance to salinity and heavy metals while maintaining key plant growth-promoting traits, including phosphate solubilization, auxin-related compound production and hydrolytic activities. In soybean (Glycine max L.), the live fungus promoted plant growth under both non-stress conditions and agriculturally relevant salinity and copper stress, and successfully colonized the roots. A cell-free extracellular extract (exCFS) promoted early seedling development in vitro, but its effects were not sustained in pot experiments. GC-MS analysis provided a qualitative characterization of the extracellular metabolite profile and revealed condition-dependent differences between control, NaCl and CuSO4; growth conditions, indicating environment-dependent modulation of fungal metabolism. Overall, these findings highlight the physiological and metabolic plasticity of T. pseudofuniculosus FRA and support its potential as a fungal biostimulant for saline- and metal-impacted agricultural soils.

Phakopsora pachyrhizi, the causal agent of soybean rust, is among the most aggressive and economically impactful pathogens affecting soybean (Glycine max) production worldwide. As a dikaryotic rust fungus, P. pachyrhizi maintains its genetic information as two discrete haplotypes that are sequestered within separate nuclei for most of its life cycle as it reproduces asexually via the production of infectious urediniospores. Previous genomes of three P. pachyrhizi isolates collected in South America revealed a large genome with a high transposable element content, measuring up to 1.25 Gb and 93%, respectively. Although these efforts provided valuable foundational resources, the complexity of the genome resulted in highly fragmented assemblies. Here, we generated a highly contiguous, substantially improved, annotated genome resource from the P. pachyrhizi isolate LA04-1, collected during the initial outbreak of soybean rust in the continental United States in 2004. The assembly is haplotype-phased, with each haplotype assembled into 18 chromosomes totaling 572 Mb and 569 Mb in size. The isolate contains unique mating-type genes and does not share haplotype genomes with any available isolates with whole-genome sequencing datasets. This resource will be valuable to future investigations into the underlying mechanisms of virulence in this important plant pathogen.

Genomic structural variants (SVs) are important sources of genetic diversity and contribute to the adaptation and evolution of plant-pathogenic fungi. Colletotrichum graminicola, the causal agent of maize (Zea mays L.) anthracnose, causes significant yield losses and serves as an important model for genetic studies. To investigate genomic plasticity in this pathogen, long-read genome assemblies were generated for four C. graminicola field strains collected in France, Croatia, Brazil, and Argentina and compared with the chromosome-level reference genome from the USA. Comparative analyses revealed extensive structural variation across core chromosomes, including large-scale translocations, inversions, duplications, and insertions and deletions, whereas minichromosomes showed markedly reduced collinearity among strains. Secreted proteins and predicted effectors were located significantly closer to regions affected by SVs. Approximately 174 kb of segmental duplications lacking genes were shared between core chromosome 04 and minichromosomes 11 and 13 in the reference genome but were absent from the other field strains, revealing a strain-specific structural link between core and minichromosomes. Protein clustering identified 733 candidate strain-specific genes across the reference genome. Approximately 72% of the genes annotated on the minichromosomes were strain-specific, indicating that most of their annotated genes lacked orthologous gene models in the other strains. Of the 733 candidate strain-specific genes, 492 retained homologous genomic sequences in at least one other strain and many also showed evidence of gene disruption. This pattern is consistent with lineage-specific pseudogenization and suggests that apparent gene presence-absence variation may arise through several evolutionary processes, including gene gain, gene loss and loss of gene function. Insertions and deletions in potentially regulatory flanking regions were associated with plant-pathogen interaction genes. In addition, interpretable machine learning analyses revealed a co-occurrence of duplications and LTR retrotransposons in genomic regions associated with secreted genes. This study shows that structural variation, transposable elements, and lineage-specific pseudogenization are important features of genome evolution in C. graminicola and identifies associations between genomic variation and genes involved in plant-pathogen interactions.

Microbial biological control agents (mBCAs) offer a promising alternative to conventional disease management, but identifying effective strains requires screening large microbial collections against many pathogens. Here, we developed a scalable screening approach to characterize the antagonistic activity of potato-associated bacteria against multiple pathogens and investigated how bacterial origin and phylogeny relate to antagonistic phenotypes. A collection of 600 bacterial strains isolated from two potato cultivars, three plant compartments using four cultivation media was screened against six major potato pathogens: the bacteria Dickeya solani and Pectobacterium carotovorum, the fungi Alternaria solani and Rhizoctonia solani, and the oomycetes Phytophthora infestans and Pythium ultimum. A scalable confrontation method was developed for non-filamentous bacteria, while complementary lower-throughput assays were used for filamentous strains. This approach enabled successful assessment of 92% of non-filamentous and 97% of filamentous strain-pathogen combinations. Isolation cultivar, plant compartment and cultivation medium had only limited effects on the proportion of antagonistic strains. In contrast, pathogen sensitivity varied markedly: P. infestans was the most susceptible pathogen, whereas D. solani and P. carotovorum were generally resistant to bacterial antagonism. Antagonistic activity clustered phylogenetically, with Bacillus and Streptomyces containing many broad-spectrum inhibitors. In addition, several taxa showed preferential activity against the oomycetes, including P. infestans-specific antagonists from the Pseudomonas genus and P. ultimum-specific antagonists from the Frigoribacterium, Curtobacterium, Pedobacter and Phyllobacterium genera. Overall, multi-pathogen screening revealed distinct generalist and specialist antagonistic profiles and identified candidate strains and taxa for further evaluation as mBCAs.

Chlorella sorokiniana is an industrially important microalga, but its genetic manipulation remains limited by poor transformation efficiency. Although several barriers have been proposed, the role of intracellular nuclease activity has received little attention. Here, we compared the degradation of exogenous double-stranded DNA in whole-cell lysates from transformation-recalcitrant (Chlorella sorokiniana and Tetradesmus obliquus) and transformation-proficient (Phaeodactylum tricornutum, Microchloropsis gaditana, and Nannochloropsis oceanica) microalgae. C. sorokiniana rapidly degraded both linear and plasmid DNA, while T. obliquus exhibited intermediate activity. In contrast, no detectable degradation was observed in the transformation-proficient species. DNA degradation in C. sorokiniana occurred within minutes, and nuclease activity in both recalcitrant species was inhibited by EDTA and Zn2+, consistent with a metal ion-dependent endonuclease. These findings suggest that intracellular nuclease activity may represent an underexplored barrier to exogenous DNA stability during transformation and provide a potential target for improving genetic transformation of recalcitrant microalgae.

Salivary secretions of herbivorous insects contain proteins that modulate plant defence and facilitate herbivory. Despite the broad host range of Spodoptera frugiperda, a major agricultural pest, its salivary molecular arsenal involved in interactions with different host plants remains poorly understood. In particular, the salivary proteins induced during feeding on different host plants and their potential interactions with plant immune components remain largely unexplored. To address this gap, we generated salivary gland transcriptomes from larvae reared on an artificial diet as well as larvae fed on tomato or chickpea plants. Comparative transcriptomic analysis identified salivary genes that were differentially expressed in response to plant feeding compared with the artificial diet. The encoded proteins were subsequently analysed using a secretory protein prediction pipeline to identify host-induced putatively secreted proteins with potential effector functions. Candidate salivary effectors were further investigated using large-scale in silico protein-protein interaction analyses with AlphaPulldown to predict their interactions with immune-pathway associated proteins of tomato and chickpea. This analysis identified multiple putative salivary effector candidates whose expression was induced by plants feeding, with distinct candidate repertoires associated with tomato and chickpea. Predicted interactions between these candidates and plant immune-associated proteins further identified potential host targets. Together, these findings suggest that S. frugiperda may adjust its salivary molecular repertoire in response to the host plant and provide a transcriptome-wide framework for validating effector-target interactions and elucidating the molecular mechanisms underlying S. frugiperda-plant interactions.

Background: Long non protein coding (lnc)RNAs are ubiquitous transcripts known to be major regulators of cellular processes in eukaryotes. In diatoms, specific lncRNA expression is tightly correlated to nutrient fluctuations, which suggests function, but their biological relevance remains to be established. In order to fill this knowledge gap, we sought to design a protocol to generate diatom lncRNA knockout mutants by CRISPR-Cas9 enabling subsequent functional studies to shed light into their mechanistic roles. Since lncRNAs do not contain open reading frames that can be targeted by point mutations for the generation of functional knockouts, we designed a multiplexed CRISPR-Cas9 plasmid harboring two independent gRNA for the specific ablation of lncRNA loci in P. tricornutum. As proof of concept, we used our plasmid-based system to target three phosphate depletion responsive lncRNAs. Results: Using the golden gate based uLoop system, we assembled a single plasmid containing all the necessary components to generate CRISPR-Cas9 large genomic ablations of lncRNAs in P. tricornutum using bacterial conjugation. Three intergenic lncRNAs (Linc91, Linc77 and Linc24) responsive to phosphate fluctuations were selected and independently targeted. Around 32% of the exconjugants presented a monoallelic ablation of the targeted loci and 11% presented a biallelic ablation on the first PCR-screening round. Mutated lines, harboring biallelic lncRNAs ablations for each lncRNA targeted, were successfully generated for each of the targeted lncRNAs. Additionally, a first line of phenotypic analysis revealed significant growth rate differences for one of the mutants ({Delta}Lnc24) under both full and phosphate depleted conditions, suggesting its involvement in cell division regulation under nutrient fluctuations. Conclusions: This work provides a simple and efficient method for the generation of lncRNA knock-out mutants in P. tricornutum. Furthermore, this approach can be easily adapted to other microalgae amenable to transformation, enabling the enhancement of our understanding of these important cellular regulatory transcripts in the marine environment.

A native entomopathogenic nematode isolate, Heterorhabditis indica Che4.1, was recovered from soils in Valencia (Spain), and its symbiotic bacterium, Photorhabdus akhurstii Che4.1, was isolated and characterized. Nematode identification was confirmed by ITS and 28S rDNA sequencing, while genome-based analyses supported the taxonomic assignment of the bacterial symbiont. The nematode-bacterium complex exhibited insecticidal activity against Spodoptera exigua and Spodoptera littoralis. H. indica Che4.1 showed an LC of 53.7 infective juveniles (IJs)/mL against S. exigua, comparable to that of the commercial reference Larvanem (54.3 IJs/mL). Against S. littoralis, Che4.1 exhibited an approximately twofold lower LC than Larvanem (77.3 versus 155.9 IJs/mL), with non-overlapping 95% confidence intervals. The complete genome of P. akhurstii Che4.1 consists of a single circular chromosome of 5.66 Mb encoding a diverse repertoire of putative virulence factors, including toxin complex (Tc) components, YenTc-like proteins, Mcf toxins, App-family proteins, Pra/Prb proteins, and chitinases, together with several predicted secondary metabolite biosynthetic gene clusters. LC-MS/MS analyses of bacterial culture supernatants and infected insect cadavers detected multiple predicted virulence-associated proteins, including Tc-related, YenTc-like, App-family, and Pra/Prb proteins. In addition, bacterial cell-free supernatants caused 68.8-81.3% mortality in S. exigua neonates under laboratory conditions. Together, these findings reveal a diverse repertoire of potential virulence determinants in P. akhurstii Che4.1 and demonstrate the insecticidal potential of its native nematodebacterium complex. The performance of H. indica Che4.1 relative to a commercial reference supports further evaluation of this isolate for biological control of Spodoptera spp., particularly under greenhouse and field conditions.

Abstract The Bacillus cereus group is a complex of closely related species, which vary in terms of their geographic spread, industrial utility, and pathogenic potential. In addition to encompassing well-known species represented by thousands of genomes (e.g., bioterrorism agent B. anthracis, foodborne pathogen B. cereus, biopesticide B. thuringiensis), the B. cereus group houses "rare" species discovered/first sequenced over a decade ago. Until recently, four rare species (B. manliponensis, B. bingmayongensis, B. gaemokensis, "B. clarus") were each represented by a single unique genome; as such, virtually nothing is known about these species beyond their type strains. We hypothesized that additional representatives of rare B. cereus group species may already be present in public databases as yet-unassembled metagenomes. To this end, we used species-specific markers to identify 7 metagenomic samples in NCBI's Sequence Read Archive, which harbored B. manliponensis, B. bingmayongensis, or B. gaemokensis. Metagenomic assembly and binning yielded an isolate-quality B. manliponensis metagenome-assembled genome (MAG) from avocado orchard soil in Israel, plus a partial B. bingmayongensis MAG from commercial spinach spiked with Shiga toxin-producing Escherichia coli. For rare species represented by >1 unique genome, we compared our novel MAGs to all available genomes and identified strain- and lineage-specific genomic elements (e.g., sphingomyelinase C and hemolysin BL homologs, variably present in B. manliponensis and B. bingmayongensis, respectively). Altogether, our study indicates that some rare B. cereus group species may be more common and geographically widespread than previously thought, and highlights the utility of culture-independent methods for recovering rare/under-isolated taxa. Impact Statement Individual B. cereus group strains vary in terms of their pathogenic potential and industrial utility. While some strains cause illness in humans and/or animals (e.g., anthrax, emetic intoxication, diarrheal toxicoinfection, non-gastrointestinal illnesses), other strains have important biotechnological applications (e.g., biocontrol, plant growth promotion, production of antimicrobials and other natural products). Because many virulence- and biotechnological-related phenotypes are strain-dependent (rather than species-dependent), phenotypes associated with an individual strain (e.g., a species type strain) cannot be extrapolated to an entire B. cereus group species. Within this context, rare B. cereus group species represent a major gap in our understanding of B. cereus group evolution, ecology, and pathogenicity, as until 2021, each rare species was represented solely by its type strain genome. Here, we aggregate all publicly available genomes/MAGs of rare B. cereus group species and supplement them with their first novel representatives since 2024. Our efforts bring the total number of unique B. manliponensis and B. bingmayongensis genomes/MAGs to >2, allowing us to conduct the first-ever phylogeographic survey and strain-level assessment of functional potential for each species. Overall, our study represents a critical first step toward understanding rare B. cereus group species beyond their type strains.

Pectobacterium comprises Gram-negative, pectinolytic plant pathogens that cause soft rot and blackleg diseases on numerous economically important crops. During sampling in 2018 in taro (Colocasia esculenta) fields in Hawaii, two bacterial strains PL152T and PL155 were isolated from taro. Peliminary molecular characterization placed both strains in the genus Pectobacterium. However, their taxanomic postions could not be resoved at the species level. A polyphasic taxonomic approach incorporating phenotypic characterization, 16S rRNA gene sequence analysis, multilocus sequence analysis (MLSA), whole genome based phylogeny, average nucleotide identity (ANI), and digital DNA-DNA hybridization (dDDH) was therefore employed to determine their taxonomic status. Phylogenetic and phylogenomic analyses consistently placed PL152 and PL155 together in a distinct and well supported lineage within the genus Pectobacterium, separate from currently recognized species. The two strains exhibited high genomic relatedness to each other, whereas the average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values, were lower than the accepted species delineation thresholds of 95% and 70%, respectively, supported the delineation of a novel bacterial species. The genomic characteristic of PL152 included a genome size of 5.02 Mb and a DNA G+C content of 52%. Collectively, based on a polyphasic approach including phylogenetic, phlyogenomic, biochemical, and phenotypic evidences supports the classification of strains PL152 and PL155 as representatives of a novel species within the genus Pectobacterium, for which the name Pectobacterium pacificum sp. nov. is proposed. Strain PL152 (ICMP 2436274, LMG 32536) is designated as the type strain.

Four rhizobial strains isolated from root nodules of tropical leguminous plants in Costa Rica were characterized. Strain LEGMi-135b was isolated from Lonchocarpus felipei root nodules, LEGMi-12c from Cojoba arborea, LEGMi-198b from Calliandra haematocephala and Rhizobium altiplani LEGMi-166a from Mimosa pigra. Phylogenetic analyses based on the 16S rRNA gene and four concatenated housekeeping genes (recA, atpD, GATase and gyrB) placed these strains within the genera Martinezella (LEGMi-12c, LEGMi-135b, and LEGMi-198b) and Rhizobium (LEGMi-166a). Genome-based analyses showed that strain LEGMi-135b exhibited an average nucleotide identity (ANI) value of 94% and a digital DNA-DNA hybridization (dDDH) value of 53.2% compared with its closest related species, Martinezella hainanensis CCBAU 57015T, supporting its classification as a distinct species. LEGMI-12c exhibited ANI and dDDH values of 98% and 87%, respectively, compared with Martinezella hainanensis CCBAU 57015T. LEGMi-198b showed the highest similarity to Rhizobium cerradonense CNPSo 3464, with ANI and dDDH values of 98% and 81.4%, respectively; however, this taxon has not been validly published. LEGMi-166a exhibited ANI and dDDH values of 100% and 99.3%, respectively, compared with Rhizobium altiplani BR 10423, confirming its assignment to this species. Additionally, differentiation of strains LEGMi-135b and LEGMi-198b from their closest phylogenetic relatives was supported by phenotypic, physiological, and biochemical analyses. Based on the combined genomic, phenotypic, and biochemical evidence, we propose the establishment of the novel rhizobial species Martinezella legmii sp. nov., with strain LEGMi-135b designated as the type strain (LEGMi-135b), and Martinezella ivanii sp. nov., with strain LEGMi-198b designated as the type strain (LEGMi-198b). This study expands our understanding of rhizobial diversity associated with nitrogen-fixing symbioses, particularly among tropical leguminous plants from Costa Rica

Background Fungi in plant root systems are key to understanding the dynamics and functioning of terrestrial ecosystems. Whereas our knowledge of below-ground plant-fungus symbioses has derived mainly from studies of mycorrhizal and pathogenic fungi, recent studies have uncovered an unexpected diversity of ecological niches occupied by root-associated fungi. The fungal order Helotiales is conspicuous in this respect because this predominantly root-associated taxon encompasses a broad ecological spectrum, ranging from saprotrophic and pathogenic fungi to endophytic and mycorrhizal fungi with plant growth-promoting functions. In this study, we integrated a worldwide dataset of 1,996 Helotiales Species Hypotheses (i.e., species-level taxa tentatively defined based on DNA sequences) detected in roots of 362 plant genera representing 117 families and 47 orders, and systematically quantified the specificity of Helotiales-plant associations. Results The analysis of 37,688 unique occurrence records of Helotiales fungi across 1,940 global sampling locations revealed highly structured associations between Helotiales genera and plant lineages, exemplified by the strong specificity of Hyaloscypha, Oidiodendron, and Phialocephala for Pinales and that of Glutinomyces for Fagales, even after statistically controlling background biogeographic patterns. At the Species-Hypothesis-level analysis, major clades within the Helotiales phylogeny consistently exhibited high specificity for Pinaceae, whereas highly specific associations with other ectomycorrhizal, arbuscular mycorrhizal, and nonmycorrhizal plants were restricted to small terminal clades of the fungal phylogeny. Conclusions The integration of ecological specificity and phylogenetic analyses suggests the possibility that Helotiales fungi established close associations with Pinaceae early in their evolutionary history and subsequently underwent host shifts to diverse angiosperm lineages. Given the emerging evidence that Helotiales fungi play pivotal ecological and physiological roles in plant root systems, these statistical insights into Helotiales-plant associations provide a fundamental basis for understanding the evolutionary history and functioning of terrestrial ecosystems.